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ABSTRACT 


Low  specific  weight  transformers  may 
be  designed  to  operate  under  impulse  conditions 
well  above  the  steady  state  limit.  This  usu¬ 
ally  results  in  nonlinear  current  and  voltage 
transformation.  A  procedure  is  outlined  in 
this  paper  for  the  analysis  of  lumped  para¬ 
meter  transformer  models  with  nonlinear  self¬ 
inductance.  An  algorithm  for  modeling  in¬ 
ductance  is  developed  which  is  accurate  for 
square  loop  cores  as  well  as  ordinary  soft 
magnetic  materials.  A  simple  routine  which 
can  be  implemented  with  modest  computing  power 
is  used  to  determine  the  dynamic  response  of 
transformers  driven  by  a  variety  of  sources. 

The  model  permits  independent  assignment  of 
initial  conditions  for  the  magnetic  state  of 
the  core  and  the  phase  of  the  driving  source. 

This  permits  a  computation  of  inrush  currents 
and  output  waveforms  under  the  entire  range  of 
possible  initial  conditions.  This  work  was 
sponsored  by  the  United  States  Air  Force  Aero 
Propulsion  Laboratory,  Wright  Patterson  Air 
Force  Base,  Ohio. 

Introduction 

The  problem  of  analysis  of  the  transient  behavior 
of  power  transformers  breaks  down  to  two  major  questions. 
First,  can  a  lumped  parameter  model  adequately  represent 
|!ie  behavior  of  such  a  device  and  second  can  a  model  of 
the  nonlinear  self-inductance  be  devised  which  has  suf¬ 
ficient  accuracy  and  flexibility  to  be  useful?  In  this 
paper,  the  validity  of  a  lumped  parameter  model  is  as¬ 
sumed  and  the  answer  to  the  second  question  is  the  subject 
of  discussion.  There  is  evidence ,  that  prior  to  the 
development  of  the  present  t  lei  the  answer  may  have  been 
negative.  Since  all  of  the  dier  circuit  elements  are 
linear,  a  Thevenin  equivale..-  for  the  network  can  always 
be  found.  Therefore  a  particularly  simple  equivalent 
circuit  is  used  to  place  emphasis  on  the  magnetic  model. 
The  resulting  circuit  equation  is  both  discontinuous 
and  nonlinear.  Thus  requiring  numerical  solution  with 
special  provisions  for  handling  the  discontinuities. 

A  simplified  lumped  parameter  equivalent  circuit 
is  shown  in  Fig.l  where  R^  ia  the  equivalent  resistance, 

Li  the  leakage  inductance,  L2  the  self -inductance  C  the 
shunt  capacitance  and  R2  the  parallel  combination  of  the 
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Simplified  trensfoneer  equivdlent  circuit 
Fig.  1 

reflected  load  resistance  and  the  effective  core  lo3s 
term.  All  of  the  components  are  linear  with  exception 

of  L2  in  this  model.  Ignoring  and  C  permits  th^ 
circuit  to  he  transformed  to  that  shown  in  Fig. 2 
where 


Vs  =  Vs 
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(2) 


R 


Typical  l-R  circi  with  nonlineer  inducUnce 
3-  i 


The  voltage  dropped  acre  the  nonlinear  inductance  L2 
is  computed  by  applying  Snell's  law  to  the  B-H  curve 
model  described  in  Ref.  1.  (  This  model  is  discussed  in 
Appendix  A.)  The  resulting  modeled  permeability  is 
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u  =  Uq 


and  the  resulting  inductive  voltage  drop  is 
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The  symbols  in  equations  (1)  through  (4)  are  identified 
in  Appendix  B.  Combining  equations  (1)  through  (4)  with 
the  remaining  circuit  voltages  in  Fig.  (2)  yields 
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This  differential  equation  could  be  solved  by  a 
number  of  standard  methods  were  it  not  for  the  dis 
continuities  which  occur  at  the  tips  of  the  hysteresis 
curves.  The  method  of  isoclines  has  been  adapted  to 
this  problem  in  a  manner  which  is  stable  under  most 
conditions.  This  technique  was  implemented  on  an 
Hp  9830  machine  using  their  modified  Basic  language . 
Returning  now  to  Fig.  (1) ,  it  may  be  shown  that  this 
circuit  is  reducible  to  the  form  of  Fig.  (2)  with  the 
result  that  the  resistance  and  source  become  complex. 
Furthermore,  the  poles  and  zeros  of  the  equivalent 
voltage  and  impedance  terms  may  be  used  to  obtain  the 
break  frequencies  for  any  given  set  of  parameter  values. 
A  specific  example  was  chosen  for  which  the  break 
frequencies  were  well  away  from  the  source  frequency. 

So  the  series  impedance  could  be  assumed  to  be  totally 
resistive.  The  parameter  values  were 


N 

=  100 

Ac  = 

4.3  E-06  m2 

£ 

=  8.47  E-  02  m 

Bs  = 

1.8  T 

Br 

=  1.0  T 

Hc  = 

40  At/m 

Fig.  (3)  is  a  plot  of  the  "ource  current  for  a  four 
volt  sinusoidal  source  wt  e  frequency  is  1  kHz  for 
R1  equal  to  0.1,  0.5,  anc.  „.0  ohms,  and  R2  =  10  ohms. 
For  this  case  as  well  as  others  in  which  a  small  core 
was  used  at  frequencies  for  which  the  other  reactances 
could  be  ignored,  the  inrush  component  in  the  current 
died  out  in  about  one  cycle.  Fig.  (4)  shows  load, 
source  and  magnetizing  currents  for  the  same  core. 
Again  the  inrush  transient  is  seen  to  disappear  after 
one  cycle. 


Concluding  Remarks 

Thi*  paper  contain*  an  outline  of  a  method  which 
is  presently  being  used  to  analyse  the  behavior  of 
steady  state  power  transformer  designs  under  transient 
conditions.  It  is  based  on  a  simple  lumped  parameter 
model  of  the  transformer  with  a  new  emperical  expression 
for  the  inductance.  The  results  of  this  technique 
applied  to  small  sample  transformers  indicates  that  for 
cases  where  the  frequency  is  low  enough  to  ignore 
leakage  reactance  and  capacitive  effects,  the  worst 
case  inrush  transient  will  die  out  in  about  one  cycle. 

The  case  for  larger  transformers  and  higher  frequencies 
have  yet  to  be  studied,  but  the  present  results  provide 
an  encouraging  indication  that  the  technique  is  successful. 
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Appei  x  A 

Magnetic  Properties  Modeling 

The  model  is  generat  from  three  experimentally 
determined  parameters.  T  ,se  are  the  coercive  force, 
the  saturation  flux  density,  and  the  residual  flux  • 
density.  The  mathematical  expressions  for  B{B)  are 
shown  in  equations  A1  and  A2  respectively. 
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K  -  J-l  (Pos.  for  uppi 
curve ) 

L  -  Inductance 
^  -  Magnetic  Path  Lei 
N  -  Number  of  Turns 
R^~  Series  Resistance 
£*2-  Shunt  Resistance 
R  -  Equivalent  Resis 
t  -  Time 

Vs -  Source  Voltage 
Vs-  Equivalent  Souro 
Voltage 


